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SUMMARY 


The objective of the study described in this report was to 
evaluate the use of coherence functions for identifying the rela- 
tive contributions of multiple dynamic inputs to the measured vi- 
bration response of spacecraft components under realistic con- 
ditions. Data for the study were generated by applying simultaneous 
vibratory and acoustic excitation of similar spectral content to 
a test model, which was a modified instrumentation truss from a Titan 
launch vehicle with a flat, lightly-loaded panel added to one bay 
to provide response data. Up to three inputs were used in tests 
conducted in an acoustically "live" room. Both independent and 
mutually coherent inputs were studied and stationary and nonsta- 
tionary environments were simulated. 

A digital computer program was written to analyze the test 
data. The mathematical approach was based on iterative computa- 
tional algorithms developed recently by Bendat. This approach is 
simpler than earlier matrix methods and is thought to be more 
economical in terms of computer time. 

It was demonstrated that the coherence function technique was 
effective in identifying and evaluating sources of excitation, for 
both correlated and noncorrelated cases, but only when the sources 
were stationary and only in the vicinity of strong structural re- 
sponses. The investigation of nonstationary sources showed that 
the application of analytical techniques based on the assumption 
of stationarity did not give acceptable results, even for rela- 
tively slowly varying inputs. 

To calculate coherence functions with high accuracy, it is 
necessary to use smoothed estimates of the associated auto and 
crosspectra. The effect of varying the number of data samples 
used in the smoothing or averaging process was briefly investigated. 


INTRODUCTION 


The situation frequently arises in the field of environmental 
dynamics in which it is desirable to identify and rank the dynamic 
sources contributing to the measured output of a system. A familiar 
case is the estimation of the relative severity of various noise 
sources making up the acoustic environment in a factory before applying 
noise reduction methods. A less familiar and more complex example 
is the problem of protecting a sensitive electronic assembly on a 
launch vehicle payload, which is being excited by simultaneous 
acoustic and vibration inputs; this can be solved more easily if 
the dynamicist knows how the various inputs contribute to the over- 
all response. In this case, the payload itself may interact with 
the inputs because of its own dynamic characteristics. 

A number of techniques are available for performing this source 
identification. If the sources are concentrated in narrow, inde- 
pendent frequency bands, frequency spectrum analysis is satisfactory 
and is frequently used. However, when the bandwidths of the various 
inputs overlap, as is often the case in launch vehicles, the re- 
sulting spectra can be ambiguous. In this situation, a time-domain 
approach can be applied, using auto- and cross-correlation functions. 
Unfortunately, there are some serious practical limitations inher- 
ent in correlation techniques, as discussed by Broch (Ref. 1); 
consequently, attention has been directed recently to using various 
forms of the coherence function to solve the problem. The coher- 
ence function is a frequency-domain function which provides the 
same basic information as the correlation function, but presents 
it in a more usable form and usually with ■ improved resolution. 

It may also be used to evaluate the fidelity of measured transfer 
functions - but this report concentrates on the application to 
solving the source identification problem. 

The purpose of the investigation described here was to develop 
a practical technique for calculating the coherence functions as- 
sociated with a structural assembly which is being excited by several 
inputs, so that the relative contributions of the inputs to the meas- 
ured vibration response could be estimated. Similar work in this 
area was reported recently (Ref. 2) in which several sources of vi- 
bration that were causing an airborne antenna to malfunction were 
investigated using coherence functions, based on a matrix treatment. 

The current study is more general and used a nonmatrix computational 
technique developed recently by Bendat (Ref. 3) which should be more 
economical than the matrix approach. Testing was performed on a 
typical structural assembly, using several input/output conditions, 
to obtain data with which to develop the technique. 
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SYMBOLS AND ABBREVIATIONS 


accel 


E { } 

dB 

f 


G. . (f) 

XX 

G (f) 

yy 

Gi.(f) 

G. (f) 
ly 

G’'f (f ) 


h(t) 

H(f) 


Hz 


X 


j 

L(f) 

MCE 

n(t) 

N(f) 

opt 

OASPL 

OCF 

P 

PCF 

Ref. 

RMS 


sec 


t 

x.(t) 


accelerometer 
expectation operator 
decibel 

frequency parameter 

acceleration unit 

auto-spectrum of x^(t) 

auto-spectrum of y(t) 

auto-spectrum of n(t) 

cross-spectrum of x, (t) with x.(t) 

1 3 

cross-spectrum of x^(t) with y(t) 
complex conjugate of G(f) 
impulse response 
Fourier transform of h(t) 

Hertz (cycles per second) 
index, i = i,2,3,...r,...p 
index, j = 1,2,3, ...r,...p 
frequency response function 
multiple coherence function 
noise input 

Fourier transform of n(t) 
optimum 

overall sound pressure level 

ordinary coherence function 

maximum value of i or j 

partial coherence function 

reference 

root mean square 

seconds 

time 

ith input 

Fourier transform of x^(t) 
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y(t) 

Y(f) 


Y 

Y 

Y 


2 

ay 

2 

y:x 

2 

iy. (i+1) (i+2) . . . 


e 


output 

Fourier transform of y(t) 

OGF between x^(t) and 
OCF between x^(t) and y(t) 

MCF between set of inputs x(t) and output y(t) 

PCF between input x.(t) and output y(t) after 
removing the effects of correlated inputs (t) 

RMS error 
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DISCUSSION OF COHERENCE FUNCTIONS 


Three forms of the coherence function are used in this report; 
these are the Ordinary Coherence Function (OCF) , the Multiple Co- 
herence function (MCF) and the Partial Coherence Function (PCF) . 

The PCF is sometimes referred to as the Conditional Coherence Func- 
tion in the literature. Other kinds of coherence functions have 
been defined (see Ref. 4) but will not be discussed here. The OCF, 
MCF, and PCF are all measures of the linear dependance between 
spectral components of random processes, and are real valued func- 
tions of frequency varying in magnitude between zero and one. 

Consider a set of stationary random inputs x^(t) (i = l,2,3,...p) 

and one output y(t) , acting on the constant parameter linear system 
shown below: 


x^(t) 

X^Ct) 


X Ct) 
p 



y(t) 


Figure 1. - Multiple input /single output system. 


The local outputs y^(t) are responses to x^(t) and are defined 
by the impulse response or weighting functions h^(x) : 


CO 

“ S x^(t-x) dx 


( 1 ) 
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The "noise” term [n(t)] shown in Figure 1 is Intended to account 
for any deviations from the ideal model; it would include measure- 
ment noise at the output, time delays in the measured data, nonlinear 
effects, statistical errors, etc. We assume n(t) to be statistically 
independent of the total measured output y(t). 

Thus , we have 


y(t) = y (t) + y (t) + + y (t) + n(t) (2) 

12 P 

Conversion of equations (1) and (2) from the time domain to the 
frequency domain by Fourier transforms results in 

Y.(f) = H.(f) X.(f) (3) 

Y(f) = Y (f) + Y (f) + + Y (f) + N(f) (4) 

1 2 P 

= H (f) X (f) + + H (f) X (f) + 


11 P P 


The spectral density functions associated with the time func- 
tions used are 




G 

yy 


(f) 


G (f) 
nn 


G. (f) 
ly 


auto-spectrum of x^(t) 
auto-spectrum of y(t) 
auto-spectrum of n(t) 
cross-spectrum of x^(t) with 
cross-spectrum of x^(t) with Y(t) 


The cross-spectrum of n(t) with y(t), G (f) , could also be de- 
ny 

fined, but would vanish since n(t) and y(t) were assumed to be sta- 
tistically independent. Note that we are using one-sided spectra, 
involving only positive frequencies. Theoretical studies often 
work with two-sided spectra, but since the correction is a factor 
of two that would always cancel out in the equations because of 
the way the spectra are used, we can safely use the more physically 
meaningful one-sided spectra immediately. 

In what follows, most quantities used will be functions of fre- 
quency, but we shall generally omit the (f) dependency for brevity 
in notation. 
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The auto-spectrum of the total output is 


yy 

But 

and 


= H 


1 1 


H 2 G + . . 

. . . + H 

2 22 

P 


H. r G, . = H. H. G. 
1 I 11 11 i: 

H.G. . = G. 

1 11 ly 

>V A 

H. G. . = G. . = G . 

1 11 ly yi 

H. 1^ G. . = H. G . 

1 11 1 yi 


so that equation (5) may be written 
P 

G = I H. G . + G 


yy 


i=l 


yi 


( 6 ) 


The ordinary coherence function between any two inputs x.(t) 
and X. (t) is given by 

^ij G. . G. . 

11 23 


(7) 


If the inputs are completely independent, the OCF will vanish 
for all combinations of i and j, since G^^ will be zero. Thus, any 

variation from zero indicates mutual coherence between the inputs. Note 

that Y?. == 1 would mean that x. and x. contain redundant infor- 
ij 1 1 

mation and one of the inputs should be eliminated from the model. 


The OCF between any input x^(t) and the total output y(t) is 



G. . G 

11 yy 


( 8 ) 


In this case, = 1 implies that the other inputs are not 

contributing to the output and the model should be viewed as a 
single input /single output system. 


If the inputs are completely independent the sum of the input- 
output OCFs should be unity for all frequencies. When this result 
is not obtained, it means that noise is present in the input or 
output data, or other inputs exist which have not been included, 
or the system is nonlinear. The OCF shows the proportion of 

the total output resulting from x^(t) ; the actual output power 

contributed by x^(t) is given directly by y? G.,. 
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The multiple coherence function shows the CF relationship be- 
tween the complete set of inputs X2(t), X2(t),... x (t) and the 

total output y(t). It can be used as a diagnostic tool, to examine 
how closely an assumed model matches the actual system. A physical 
interpretation of the MCF concept may be obtained as follows: 


Let y(t) = predicted output from an ideal model 
y(t) = actual measured output 

A. 

n(t) = effects causing deviations = y(t)-y(t) 

G and G = corresponding auto-spectra, so that 
yy yy nn t- & t' » 

G = + G 

yy yy nn 


( 9 ) 


The MCF is defined as 
G^^ G -G 

^2 = _ZZ = nn 

^x:y G G 

yy yy 


( 10 ) 


Since both G"^ and G are real and positive-valued, G can 
yy nn ^ yy 

2 

never be less than G^^ so that the inequality 0 < Yy.^^ < 1 is always 

satisfied. Under ideal conditions, the MCF will be identically 
equal to unity for all frequencies. 


If all the inputs are independent, the MCF will equal the sum of 
the individual input-output OCFs: 





+ . 



( 11 ) 


This equation does not apply if the Inputs have some degree of 
mutual coherence; however, a similar relationship can be developed 
between the MCF and various partial coherence functions . Before 
the PCFs can be calculated the inputs must be modified or "condi- 
tioned" to remove the coupling effects caused by the lack of in- 
dependence between the inputs. Each PCF will then show the pro- 
portion of the output caused by a specific input, with the effects 
of all other inputs eliminated. 


Special notation is used to show that the data records have 
been conditioned. Thus, for example, X3.22(t) indicates that in- 
put X3(t) has been conditioned on x^(t) and X2(t) to eliminate 
their effects. The corresponding auto-spectrum would be G33 22* 

The auto-spectrum of the output after conditioning on X2(t) and 
Xo(t) would be G . The PCF between Xo(t) and y(t) , conditioned 

yy-12 

on X2 (t) and X2(t), is given by 
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2 


Y 


2 

3Y. 12 


= 1^37. J- 2 I 

-XT C 

33.12 yy ■ 1 2 


Thus, the PCF between the preconditioned quantities is equiv- 
alent to the OCF between the postconditioned quantities. 


The numerical conditioning process is quite complicated for a 
multiple input system. Matrix methods have been developed (Refs. 

4, 5, 6) and an application is described in reference 2. These 
methods are effective in keeping the computational process orga- 
nized, but can be expensive in computer time, involving a consid- 
erable number of complex matrix inversions. An alternative method 
was developed recently by Bendat(Refs. 3, 7), in which iterative 
computational algorithms are used to calculate conditioned spectral 
density functions and partial coherence functions. The technique 
uses Fourier transforms of data records and performs the condition- 
ing by removing an optimum least squares prediction of the "contam- 
inating" inputs from each input in turn. A unique feature of 
Bendat’s method is the use of "ordered" inputs. In this process, 
the original set of inputs records are rearranged in a particular 
order and then each record is conditioned on the preceding records. 
Thus, Xi(t) is unchanged, X£(t) is conditioned on x^Ct) to become 
^2.lX*^)> is conditioned on x^Ct) and X 2 (t) to become 

and so on. The noise term n(t) is now equivalent to the output 
conditioned on all of the inputs. The output term y(t) remains 
in the model in its original nonconditloned form, however. 


A revised model involving ordered, conditioned, frequency-domain 
inputs is shown in Figure 2. 



Figure 2. — Multiple input /single output system with 
ordered conditioned inputs. 
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The following procedure is used to establish the new order 
for the inputs: 

(1) Calculate the OGF between each input xj(t) and the output 
y(t). 

(2) For the frequency range of primary interest, arrange the 
calculated OCFs in order of magnitude. 

(3) The new xi(t) is the input which gives the highest OCF, 
X 2 (t) is the input giving the next highest OCF, and so on for the 
remaining inputs . 

This gives a set of ordered inputs which now have to be condi- 
tioned by the process developed in reference 3 and summarized here. 


The block labeled (f ) in the model represents frequency re- 
sponse functions connecting the conditioned inputs x^ ^ (f ) 

with each other and with the output Y(f). The are used to 

estimate the effects of the preceding (i-1) inputs which are to be 
eliminated from the i th input in the interactive conditioning proc- 
ess. The L,, must be optimized on the basis of minimizing the 
ij 

error in these estimates. 


From the model, we can write 
P 

Y = X! L. X. + N 

ly 1.123 (i-l) 

Let N. = Y- L. X. 

ly ly 1. 123* • • 


( 12 ) 

(13) 


which is the difference between the total output y and the output 
due to the i th conditioned input passing through 1‘^y* 


Then 


N 


ly 


^ = |y - L. X. 


A 


iy ^i* 123* • • (1-1) 


= Iy 2 - L. [X. f. T . Y 

I I ly 1. 123 • . ♦ • (i“l) ] 

* , A 

-L. [YX. L. 

ly 1. 123- . . (i-l) ly ly 


123 


..(i-1) 


(14) 
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- where the asterisks indicate complex conjugates. 

Taking expected values of both sides gives an expression for the 
mean square system error for any terms of auto-and cross- 

spectra of the conditioned input and the output: 


E{IN. |2} = e2 
' ly ‘ 

= G -L. '''g. .... 

yy ly ly. 123 * • • (i“D 

A A 

^iy ^ iy. 123 • • • (i"l) ^iy ^iy ^ii. 123 • • • (i“l) (15) 

To minimize e^, differentiate equation (15) partially with re- 
spect to L, and equate to zero: 


i(e^) 

a-L. 

ly 


iy . 123 • • 


V-L-Xy 




Taking complex conjugates and solving for the optimum gives 


L^y (opt) = iy. • • ( 1 - 1 ) 

‘^ii- 123... (1-1) 


( 16 ) 


The optimum input/output frequency response functions for the 
model will thus be 


Ll 


L 2 


^ Giy 

y — - 

Gi 1 
^2y • 1 


y G 

22. 1 


l3 


y 


Gsy • 12 
G33. 12 


etc . 


In a similar manner, we can find the optimum frequency response 
functions connecting the coherent inputs, by setting y=j : 


L_ (opt) 


- ij . 1 9 . (i- 1 ) 

^11 • 123 • ' • (i“l) 


i = 1 , 2 , p 

j > 1 


( 17 ) 
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I 1— J 
1 

Li o = 

3 

1 1 

Gi 1 

23. 1 

L = ^24.1 

^24 


^ 22.1 ^ 22.1 


Q 

^(p-l)p " — (-£-~A)p_- ■ V (P - 2) 

^(p-1) (p-1) . 123- • • Cp-2) 

At this point, we have developed expressions for the optimum 
frequency response functions in terms of conditioned- auto- and 
cross-spectra; now we need to be able to obtain these spectra from 
the original nonconditioned spectra. An iterative algorithm was 
derived by Bendat to do this. (The derivation is lengthy and will 
be omitted here. Details may be found in reference 3.) The gener- 
al expression to be used is: 


G. . 

IJ • 123 


r 


G 


ij . 123 . . . (r-1) 


L .G. 

rj ir.i23 


. . (r-1) 


where i=l , 2 , 3 , . . . p ,y ; j £ i; r < j 
Special cases are 


r=l: 


r=2: 


r=3: 


G.. =G.,-1. ,G. 

iJ-l iJ IJ 11 


i,j 


G.. =G.. -L.G, 1 

iJ-12 iJ-1 2l 12* 

G.. =G.. -L.G. 

lJ-123 1J‘12 3J 13- 


= 2,3, 
i,j = 


• -p>y 


3,4, . . .p,y 
= 4,5,6, . . .p,y 


(18) 


r=p: G 

yy • 1 2 3 • • • P 


G 


yy • 123 • • • 


L 

py 


G 


yp. 123- . . (p-1) 


To obtain the auto-spectrum form, set i=j in equation (13) : 


G. . 

33 ■ 123 


^ •123‘ • • (^-1) 


L.G. . 

rj jr. 123- • • (r-1) 


Now rewrite equation (17) with i replaced by r: 


L . 

rj 


G . 

= rj 


123 


....(r-1) 


G 

123 


..(r-1) 


(19) 


(20) 
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, then 


Take complex conjugates and solve for G / 

jr . 123 • • • (r-1) 


substitute in equation (19) to obtain 


G. . 

JJ -123- • 


'j j • 123- • • 


2 

^rjl ^rr. 123- . . (r-1) 


( 21 ) 


Again, we could evaluate special cases by setting r=l,2,...p and then 
varying j so that j=r+l, r+2,...p,y. The iterative process of obtaining 
the ordered conditioned auto-spectra from the original auto-spectra is 
illustrated in Figure 3. 


11 


22 


33 


"44 


55 



S 5 .I ^55.12 S 5 .I 23 S 5 .I 234 


Figure 3. - Diagram showing how ordered conditioned atuo-spectra are 
obtained from original auto-spectra (from Ref. 3). 


From equation (21) it can be proved that, for any r < j where r j< p. 


r 

G. . = > |L. . G. . . . + G. . 

33 X— rfl xj I XX. 123- • • (i-l) JJ-123*-*^ 

i=l 


( 22 ) 


Now if we set r=p and let j=y, equation (22) yields a form similar 
to equation (6) : 


r 

^jj "S l^iyT ^ii.i23--*(i-l) ^nn 

i=l 


(23) 
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The PCF between the conditioned inputs and outputs are now 
defined as 

|G. I ^ 

^2 I ly. . . (x-1) I 

iy • 1 2 3 • • • ( G 

xi.123- • • (i-l) yy- 123- • • (l-l) 


( 24 ) 


A general formula connecting multiple coherence functions with 
associated ordinary and partial coherence functions is derived in 
reference 7. For a system with p inputs, this formula is 


^ ^--^2y.P *^^-^3y.l2^ ' • ' • ^^-Ypy.123. . . (p-1) 


)} 


(25) 


Note that = 1 if and only if one of the OCF and PCFs equals 

y:x ^ 

unity, and = 0 if and only if all of these functions equal zero. 

y;x — 

For the case where the inputs are independent, equation (25) reduces 
to the simpler result in equation (11) . 
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ANALYSIS OF THREE INPUT/SINGLE OUTPUT SYSTEM 


The general theory discussed in the previous section was applied 
to the special case of a system with three inputs, which may either 
be independent or mutually coherent, and a single output. This sys- 
tem corresponds to the test setup used to obtain data from which the 
computational technique was developed. 


The test results provided inputs Xi(t), X2(t) , X3(t) and an out- 
put y(t) ; both independent and coherent inputs were applied as sepa- 
rate cases. These time histories were then processed to yield averaged 
estimates of auto-spectra Gn(f), G22(^)» ^33(1) and G (f) and (for 


the coherent case) cross-spectra Gi2(f)> Gi3(f), Gj^^(f) , G23(f), G2y(f) 
and G3^(f). For the independent input case, the OCFs were computed 


using equations ( 7 ) and (8), and then the MCF calculated from equation 
( 11 ) . This was the extent of the computations performed for the in- 
dependent input case. In the coherent input case, the OCFs were com- 
puted for use in reordering the inputs. The ordered, conditioned 
auto- and cross-spectra were next calculated from equations ( 18 ) and 
( 19 ). To start the iterative computation process, it is necessary to 
calculate the first set of optimized frequency response functions from 
equation ( 17 ) : 


^12 





( 26 ) 


- note that these are ratios of nonconditioned spectra only. 

The relationships in equations ( 18 ) , ( 21 ) , and ( 26 ) allow us to 
calculate the following: 


G22. 1 

- G22 “ 

1^12!^ Gii j 


^ 33 . 1 

G33 - 

1 ^13 ! ^ 1 j 

( 27 ) 

yy.l 

= G 

yy 

1 ^ly|^ ^11 1 


G23 . 1 

* 

= G23 “ Li 3 G 12 1 


G2y. 1 

^ I 

= G2y - Li G 12 1 

1 

( 28 ) 

^ 3 y. 1 

1 

•k 

= Gsy Liy G 13 
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The next set of frequency response functions can now be calculated: 

L23 = Lay = (29) 

^ 22.1 ^ ® 22.1 

Knowing these, we can compute 


G 33 . 

. 12 

“ G 33 _ ] 

1 ^ 23 ] 

i 2 p 

G 22.1 


G 

yy< 

, 12 

" G 2 y . 1 

" l^zyl 

^ ^ 22. 1 1 

1 

! (30) 

G3y- 

' 12 

= Gsy.i 

* 

- L 2 y G 23 . 1 


Next , 






Lsy 

G 

3y. 12 



(31) 

G 

33*12 





Finally, the auto-spectrum of the noise at the output may be calculated: 


G = G 

yy -123 nn 

^ Gyy “ |^ly|^ ^11 " |L2y|^ G22-1 - |^ 3 yj ^ G33.12 ( 32 ) 


Using equations (27), (28), and (30), 
tions may be calculated: 

y 2 _ |G2y.i|^ 


2y.l 


G 22 . 1 


G 

yy - 12 


the two partial coherence func- 


(33) 


3y. 12 


r ^ 

'^3y . 12 
^33 . 12 


(3A) 


yy. 12 

The multiple coherence function is given by 
^y:x" ^^~^2y.l^ ^^"^3y.l2^ 


(35) 


The computational process is summarized in the flow chart, 
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EXPERIMENTAL PROGRAM 


Testing was performed on a typical aerospace structure to 
provide data for use in demonstrating the computational technique. 
The structure was an instrument support truss from a Titan launch 
vehicle, shown in Figure 5. A panel, carrying a simulated compo- 
nent, was mounted on the top of the truss (see Fig. 6). It was 
designed so that its first three resonances would occur in the 
range 50 to 200 Hz. The assembly was supported inside the Martin 
Marietta acoustic chamber on a 2 Hz suspension system and was 
connected via two push rods to two electromagnetic shakers. The 
shakers were mounted on rubber isolation pads to minimize feed- 
back of mechanical vibration through the floor. 

The test setup was designed to allow application of two sep- 
arate random mechanical vibration inputs through the truss to the 
panel, combined with direct acoustically-induced vibration. Thus, 
as many as three inputs could be applied simultaneously. The out- 
put was defined as the response measured at the simulated component 
on the panel. 

The test program was conducted in two phases. In Phase I, all 
of the inputs were stationary; that is, their statistical properties 
were time invariant . All of the theoretical developments discussed 
earlier are based on the assumption of stationarity . In Phase II, 
the inputs were deliberately made nonstationary so that, strictly 
speaking, the definitions of the various coherence functions were 
no longer valid. The data provided by Phase II were required so 
that the errors introduced by falsely assuming stationarity could 
be evaluated for a particular case. 

Table 1 lists the test runs that were performed in Phase I 
showing the various input conditions. For the "coherent input" 
case, the tv70 shakers were driven through a single power amplifier/ 
random noise generator system. For the case of three independent 
inputs, separate power amplifiers and random noise generators were 
used . 

Table 2 shows the inputs used for the Phase II test. Only two 
inputs were used, consisting of mechanical vibration from one shaker 
combined with acoustics. It was considered that this setup was 
sufficiently complicated, in view of the additional complexity 
introduced by the lack of stationarity. 




Figure 5. - Test setup 









TABLE I. - PHASE I TEST RUNS 



Input Conditions 


Run No . 

Shaker 1 

Shaker 2 

Acoustxcs 

Notes 

1 

X 

X 


Shakers coherent 

2 

X 

X 

X 

Shakers coherent, 
acoustics indept . 

3 

X 

X 

X 

Indept. inputs 

4 

X 




5 



X 


6 

X 


X 

Indept. inputs 


TABLE 2. - PHASE II TEST RUN 



The transducers used to measure the inputs and the outputs are 
identified in the notation used for the data analysis, in Figure 7. 
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DISCUSSION OF RESULTS 


Phase I: Stationary Inputs 

Checks were first performed to verify that the digitized data 
were correctly scaled. The time history from one of the acceler- 
ometer channels was played back on an oscillograph record direct- 
ly from the analog tape and compared with the corresponding digi- 
tized time history. Next, the acceleration spectral density was 
plotted using an analog data analysis system and compared with 
the acceleration spectral density obtained from the digitized time 
history using a Fast Fourier Transform (FFT) method. Finally, a 
one-third octave band Sound Pressure Level spectrum was plotted 
from a sample of the analog microphone data and used to calculate 
the corresponding Pressure Spectral Density plot; this was then 
compared with the same plot resulting from the digital analysis. 
The three comparisons confirmed that the scaling used in the ana- 
log-digital conversion and the subroutines written to calculate 
the spectra were correct. 


The technique used on the test data consisted of first calcu- 
lating the auto- and cross-spectra for the inputs and output sig- 
nals, averaged over 50 segments of data, then computing the OCF 
between each pair of inputs (y^-) ^.nd between each input-output 

combination (y^ ) . The y^ . were then reviewed to determine whether 

ly 13 

the two inputs and x, should be treated as being independent 

or correlated (i.e., mutually coherent). According to Bendat (Ref. 
8), a value of y^^ of 0.1 is low enough for an assumption of inde- 
pendence; a higher value indicates significant coherence between 

the inputs. At the other end of the scale, y^ . ^ 0.9 means that 

ij 

the two inputs are effectively fully coherent and should therefore 
be treated as a single input to the system , 


For the runs in which the inputs were judged to be indepen- 
dent, the MCF between the total set of inputs and the output was 
then computed by adding up the y^ , as indicated by equation (11) . 

For the tests involving coherent inputs, the y?^ were used to 

reorder the inputs (if necessary) as discussed earlier and then 
PGFs and the MCF were computed, using equations C26) through (35). 
The results of treating the data by these techniques will now be 
discussed for each test run. 
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Run 1 (two shakers, correlated ) . - Acceleration spectral den- 
sity plots for the two inputs and the response are presented in 
Figure 8. The inputs are similar in spectrum shape at both shakers, 
and reflect the resonances and antiresonances of the truss struc- 
ture. The response of the panel is shown to occur as prominent 
peaks at approximately 68 Hz, 161 Hz, 230 Hz and 410 Hz, with many 
intermidiate minor peaks. 

Coherence function plots are given in Figures 9 through 11. 
Figure 9 is the OCF between the two inputs. It shows that the 
two inputs were well correlated; in fact, by Bendat's criterion, 
their mutual coherence was high enough (>0.90) over much of the 
frequency range that the system should be regarded as having only 
a single input. On this basis, the OCF between either input and 
the response should be essentially the same, and should theoreti- 
cally be unity across the frequency band if the system were linear 
and noise-free. A comparison of Figures 10 and 11 shows that they 
are very similar; however, neither OCF ever reaches unity, indi- 
cating the existence of nonlinearities and/or measurement noise 
in the system. 

Run 2 (two shakers, correlat ed, plus acoustics) . - Accelera- 
tion and pressure spectral density plots for the three inputs 
and the response are given in Figure 12. The shaker inputs were 
again very similar in spectrum shape and the response peaks oc- 
curred at approximately the same frequencies as in Run 1. 

The OCF between the two shaker inputs is plotted in Figure 
13. This shows that the coherence between the inputs was consid- 
erably less than for Run 1, particularly in the frequency range 
below 100 Hz. However, the OCF was in the desired range of 0.10 
to 0.90 over most of the frequency band. The OCFs between each 
individual shaker input and the acoustics input were examined next; 
these are presented in Figures 14 and 15. Although the shaker in- 
puts and the acoustic inputs were theoretically independent, a 
significant degree of mutual coherence is indicated. This is 
thought to be caused by the acoustically induced panel response 
being transmitted through the truss structure and coupled into 
the vibration environment as recorded by the shaker input acceler- 
ometers. The Run 2 data were therefore treated as representing 
a system with three correlated inputs and a single output. 

The resulting coherence functions are shown in Figures 16 
through 20. When these are examined in conjunction with the panel 
response plot [Figure 12(d)], several features become apparent: 

(1) Figure 16 indicates y^. = 0.78 at 68 Hz, showing that 
the high response at this frequency resulted mainly from the shaker 
1 input. Figures 17 and 18 are plots of the Partial Coherence 
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ACCEIZRATION SPECTRAL DENSITY - g'/Hm 
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Cb) Shaker 2 input 



Figure 8: Input and response data for run 1 
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Ordinary coherence function 



Frequency - Hz 


Figure 9: Ordinary coherence function between 

shaker #1 and shaker #2 inputs (run 1) 
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Ordinary coherence function 
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Figure 10. - Ordinary coherence function between 
shaker 1 input and panel response 
(run 1) 
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Ordinary coherence function 
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Figure 11. - Ordinary coherence function between 
shaker 2 input and panel response 
(run 1 ) 
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Figure 12. - Input and response data for run 2 
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"Figure 14. - Ordinary coherence function between 
shaker 1 and acoustics (run 2) 
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Figure 15 - Ordinary coherence function between 
shaker #2 and acoustics (Run 2) 



Ordinary coherence function 



Figure 16. - Ordinary coherence function between 
shaker 1 input and panel response 
(run 2) . 
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Partial coherence function 
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Figure 19. - Multiple coherence function between total 
set of inputs and panel response (run 2) . 




Functions (PCF) between shaker 2 and panel response with the ef- 
fects of shaker 1 removed, and between the acoustic input and 
panel response with the effects of both shakers removed. Both 
PCFs have low values at 68 Hz, indicating little contribution to 
the response. 


(2) The next response peak is at approximately 160 Hz. Fig- 
ure 16, 17, and 18 show that the two shakers were primarily re- 
sponsible for this peak, having OCF and PCF values of 0.44 and 
0.47 respectively. The PCF associated with the acoustic input at 
this frequency is less than 0.10. It is interesting to note that, 
from Figure 13, the mutual coherence between the two shakers was 
very high (aj; 0.98) in the 150-180 Hz. range, whereas the mutual 
coherences between the Individual shakers and the acoustic input 
(Figures 14 and 15) were close to zero (^0.02). Thus, in this 
frequency band, the system should actually be analyzed as a model 
having two independent inputs and a single output . 


(3) The third response peak, occurring at about 230 Hz, was 


due mainly to the input from shaker 2, since y 


Y 


2 _ 


2y.l 


= 0.53, whereas 


= 0.15 and ^2 ~ this frequency. 


(4) Another narrow response peak is seen at approximately 

290 Hz. The CF values at 290 Hz are y^ = 0.02, Yo . = 0.01 and 

ly 2y.l 

Y^^ ~ 0.85; thus, this peak was clearly caused by the acoustic 

input . 


(5) A rather broad response peak is centered at about 420 

Hz. This also appears to have been caused by the acoustic input, 

which has y^ = 0.6. The OCF and PCF for the shakers are ap- 

3y.l2 

proximately Y^y = 0.05 and Y^y “ 0.06. 

The following table3 summarizes the ordinary and partial coher- 
ence functions associated with the first six response peaks on the 
panel. The multiple coherence function values (taken from Fig. 

19) at these frequencies are also tabulated and compared with the 
sum of the three individual coherence functions. 


The sum of the individual coherence functions is close to unity 
in every case, and is higher than the MCF at all of the modal fre- 
quencies. It can easily be shown that will always exceed 

Y^ , as follows: 
y :x 

From equation (35), Y^.^ = 1 - (1 - Y^y) (1 ~ ^2y.l^^^ ~ *^3y.l2^ 
for a three input system. For brevity, write this as 














( 36 ) 


Yy.^ = 1 - (1 - A) (1 - B) (1 - C) 

=1- (1-A-B-C+AB+BC+CA- ABC) 

= (A + B + C) - (AB + BC + CA) + ABC 

In the same notation, + B + C (37) 

= (AB + BC + CA) - ABC 

= AB (1 - C) + BC + CA (38) 

Since | A, B, C } ^ 1.0, the quantity (1 - C) is always positive 
or zero, so that Ey" must always be greater than the multiple co- 
herence function 

y:x 

The amount by which the MCF is less than one represents the 
effects of noise and nonlinearities in the system. Since these 
effect's may be assumed to be statistically independent of the in- 
puts they do not contaminate the individual coherence functions 
associated with the inputs, only the multiple coherence function, 
which represents the coherence between the full set of inputs and 
the measured output. 


As a final comment, Figure 20 allows a direct comparison of 


the PCF between the acoustic input and the panel response (y 


3y.l2, 


Fig. 18) with the OCF that is obtained when the mutual coherence 
between the acoustic input and the shaker inputs is neglected. 

It can be seen that the difference between the two CFs is very 
small, implying that the errors introduced by the simplified ap- 
proach would be negligible in this particular case. 


Run 3 (two shakers, noncorrelated , plus acoustics) . - The in- 
put and response data are plotted in Figure 21. It was intended 
to use shaker inputs that were uncorrelated with each other as 
well as with the acoustic input, giving three independent inputs. 
However, when the OCF between the two shakers was plotted (Fig. 
22) , it was found that high mutual coherence existed at some fre- 
quencies between the two measured shaker inputs. This was es- 
pecially true in the 130-180 Hz region, where a lateral mode ex- 
ists for the truss, allowing the two shaker input locations to 
be very closely coupled. To a lesser extent, the acoustic input 
was found to be coherent with each individual shaker input (Fig- 
ures 23 and 24) ; therefore, it was decided to treat the run 3 
data as applying to a system with three correlated inputs and a 
single output, like the run 2 data. 


The coherence functions between the three inputs and the out- 
put are given in Figures 25 through 28 . The plots can be used 
to identify the contributions of the three excitation sources at 
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Ordinary Coherence Function 



Frequency - Hz 


Figure 22. - Ordinary coherence function between 

shaker 1 and shaker 2 inputs (run 3) . 
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Ordinary coherence function 



Figure 23. - Ordinary coherence function between 
shaker 1 and acoustics (run 3) . 
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Ordinary coherence function 




Ordinary coherence function 


0.70 



Figure 25. ~ Ordinary coherence function between 
shaker 1 input and panel response 
(run 3) . 
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Partial coherence function 
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Figure 26. " Partial coherence function between 
shaker 2 input and panel response 
(run 3) . 
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Figure 27- - Partial coherence function between 
acoustic input and panel response 
(run 3) . 


Multiple coherence function 
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Figure 28. - Multiple coherence function between 

total set of inputs and panel response 
(run 3) . 
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the various panel resonances, in the same way as the CFs for the 
run 6 data were used. Table 4 summarizes the values of the ordi- 
nary, partial and multiple coherence functions which were taken 
from the plots. 


TABLE 4. - VALUES OF COHERENCE FUNCTIONS AT PANEL RESONANCES (RUN 3) 


Frequency, Hz 

ly 

2y.l 

Y ^ 

3y.l2 

1Y2 

Y2 

y:x 

68 

0.37 

0.39 

0.37 

1.13 

0.77 

161 

0.58 

0.26 

0.01 

0.84 

0.69 

415 

0.27 

0.26 

0.40 

0.93 

0.67 

580 

0.23 

0.50 

0.46 

1.19 

0.79 

810 

0.11 

0.37 

0.63 

1.06 

0.82 


Note that the ZY^ values are again close to unity, and in 
three cases actually exceed unity. This probably indicates that 
errors are affecting the individual CF values contributing to 
the summation. For example, at 68 Hz, th.e inputs from Shaker 1 
and the acoustics are essentially independent, since Yin is only 
about 0.05. Also, at 161 Hz the CF between the two shakers is 
greater than 0.9, so that they should probably be treated as a 
single input in this frequency band. The conclusion is that, to 
obtain the most accurate results, several different interpretations 
of the system would be necessary, changing the model as appropri- 
ate for different frequency ranges as indicates by the mutual co- 
herence between the inputs. 

R un 4 (shaker 1 only) . - The objective of this run was to sup- 
plement the data from Runs 2 and 3 and to identify the amount of 
"noise” in the system. The OCF for the single input/single out- 
put system would be 1.0 across the full frequency range in the 
absence of noise or nonlinearities. Figure 29 shows the input and 
response data. The plots are very similar to those obtained from 
Run 4 which, as discussed earlier, was effective as a single in- 
put/single output case. The response data plot shows high peaks 
at 68 and 161 Hz. These were the frequencies at which the shakers 
were responsible for most of the response, according to the in- 
terpretation of the associated PCFs. Thus, the Run 4 data sup- 
port that interpretation. Figure 30 is a plot of the OCF and is 
similar to the OCF plots (Figures 10 and 11) from Run 1. These 
show that the analyzed data have a high noise level at most fre- 
quencies in the range 5 to 1000 Hz. Low damping and, hence, low 
off-resonance response and narrow frequency bands contribute to 
this problem. 
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Acceleration spectral density - g^/Hz 




(b) Panel response 


Figure 29. - Input and response data for run 4 . 
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Ordinary coherence function 




The Run 4 data were also used to investigate the effect of 
varying the number of data segments used in the averaging process 
before calculating the OCF. Figure 31 shows four versions of the 
OCF, in which 10, 50, 100, and 200 averages were used. With only 
four different cases, the results were not conclusive; however, 
as illustrated in Table 5, the variation of some frequencies was 
large enough to suggest that a number of averages in the region 
of 200 is necessary to ensure reasonable accuracy. On the other 
hand, if a fixed record length of data is available for analysis, 
the accuracy of the results will be degraded if the length of 
each segment is reduced, so a detailed error analysis is necessary 
to ensure that the best trade-off is made. 


TABLE 5. - EFFECT OF AVERAGING ON OCF VALUES 



Finally, data from Run 4 were used to plot the OCF in a 
slightly different format, using linear scales for both coherence 
function and frequency. The result is shown in Figure 32 and is 
based on using 200 averages. The frequency resolution is worse 
at low frequencies and better at high frequencies than the cora-- 
parable semilog plot, and the preferred format will depend on the 
particular application. 
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Figure 32. - Ordinary coherence function between shaker 1 
input and panel response, using 200 averages, 
with linear frequency scale (run 4) . 
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Run 5 (acoustics only) . - This run, like Run 4, was intended 
to supplement Runs 2 and 3. Since the panel response in this run 
was caused only by acoustics, it can be compared with the acous- 
tically-induced responses in Runs 2 and 3, as indicated by the 

PCFs CY^ 1 o) • the discussion of the Run 2 results, it was 

3y.l2 

concluded that the responses at 290 and 420 Hz were caused pri- 
marily by the acoustics. The response plot in Figure 33(b) sup- 
ports this conclusion strongly at 420 Hz, and to a lesser extent 
at 290 Hz. The OCF plot in Figure 34 again shows that significant 
noise and/or nonlinearity effects contaminated the measured data. 


Run 6 (shaker 1 plus acoustics) . - Two nominally uncorrelated 
inputs were used for this test run. Figure 35 shows the input 
and response data. As usual, the coherence between the inputs 
was first calculated to verify that the assumption of independence 
was correct. Figure 36 shows that this was the case except at 
the 4.8 Hz point and at some frequencies above 500 Hz when the 
value of the OCF exceeded the Bendat criterion of 0-10 by a small 
margin. It was concluded that it would be reasonable to ignore 
this degree of mutual coherence and treat the system inputs as 
uncorrelated . 


The OCF plots between each input and the output appear in Fig- 
ures 37 and 38. When these are reviewed in conjunction with panel 
response plot in Figure 35(c), it is clear that the response peaks 
at 68 and 160 Hz are due primarily to the shaker input, while the 
resonance at 420 Hz appears to be driven most efficiently by the 
acoustics. The MCF plot in Figure 39 again shows that a high de- 
gree of "noise" is present in the data. 
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(a) Acoustic input 


(b) Panel response 


Figure 33. - Input and response data (run 5) 
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Ordinary coherence function 
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Figure 34. - Ordinary coherence function between 
acoustic input and panel response 
(run 5) . 
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(a) Shaker 1 input 



(b) Acoustic input 



Figure 35. 


Input and response data (run 6) . 
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Figure 36. - Ordinary coherence function between 

shaker 1 input and acoustics (run 6) . 




Ordinary coherence function 


F 



60 


Ordinary coherence function 



Multiple coherence function 





Phase II: Nonstationary Inputs 


Data from Run 7 were used to investigate the application of 
the CF approach to a system with two obviously nonstationary un- 
correlated inputs. The inputs were a single shaker plus acous- 
tics, and were initially identical to those used in Run 11 (see 
Fig. 35) . After applying the stationary inputs for a few seconds, 
the overall acceleration and sound pressure levels were gradually 
increased to peak values approximately 8 dB above the initial 
values, and then reduced back to the initial values. The time 
variation of the overall levels was as shown below: 


Overall 

level 



To obtain a basis for comparison, the OCFs between the indi- 
vidual inputs and the response were calculated using data from the 
first 14 seconds, when the inputs were stationary. 

The roll-off part of the time history consisted of the last 
30 seconds of the test, equivalent to 150 data segments after dig- 
itizing. This was used for the nonstationary analysis. Six dif- 
ferent data records, each containing 20 data segments, were ex- 
tracted by the following method: 

(1) For the first data record, the first 20 segments out of 
the 140 available were used to calculate two OCFs. 

(2) For the second record, alternate segments out of the 
first 40 were employed for the calculation. 

(3) The third record was made up by taking every third seg- 
ment until 20 segments were accumulated. 
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(4) For record #4, every fourth segment from the first 80 
segments was used. 

(5) The fifth and sixth records were made up in the same way, 
utilizing every fifth and sixth segment respectively. 

The six different analyses thus covered a range in record 
lengths of approximately 4 to 24 seconds, with a corresponding 
spread of about 6 dB (factor of 2) on the overall levels of the 
inputs . 

Figure 40 shows the OCF between the shaker input and the panel 
response for the stationary precursor in Run 7. Figures 41 and 
42 give the same OCF for the six increasingly nonstationary time 
records. The OCF between the acoustic input and the panel re- 
sponse for the stationary case and the six nonstationary cases 
is plotted in Figures 43, 44, and 45. 

A comparison of for the stationary part of the run, plot- 

ted in Figure 40, witn’^he same OCF computed from the six nonsta- 
tionary time records (Figures 41 and 42) shows that the general 
shape of the OCF is similar for all seven conditions. However, 
large variations occur in the magnitude of the OCF, with no ap- 
parent relationship to the length of the time record. An inspec- 
tion of the plots for y^ leads to the same conclusion; that the 

3,y 

introduction of nons^ationarity to the inputs causes tne OCR to 
change significantly, but not in a predictable way. 
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(d) Data record 4 


(e) Data record 5 


(f) Data record 6 


Figure 42. - Ordinary coherence function between shaker 1 input and panel 
response, nonstationary part of run 7 . 
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CONCLUSIONS 


Phase I: Stationary Inputs 

From the results of Phase I of the study, it was concluded 
that the coherence function technique is effective in identify- 
ing the relative contributions of multiple sources of dynamic ex- 
citation to the response of a practical mechanical system. How- 
ever, this was found to be true, for the system tested, only in 
the region of resonant response frequencies. The technique did 
not give good results at frequencies between resonances, when er- 
rors induced by noise and/or system nonlinearities were apparent- 
ly similar in magnitude to the quantities being used in the cal- 
culation. This is not a serious handicap when using the techni- 
que strictly for source identification, since such information is 
generally of greatest interest at resonant frequencies, but may 
seriously hinder applications in which general inputs are being 
deduced from the response of a specific resonant system. The 
noise levels were high in all of the test cases, as indicated by 
the fact that the multiple coherence function approached unity 
only at a relatively few frequency points. For an ideally linear 
and noise-free system, the MCF should have a value of 1.0 across 
the full frequency band if all the appropriate inputs have been 
properly included. 

The iterative approach developed by Bendat (Refs. 3 and 7) 
was used to calculate the partial coherence functions for the 
cases involving correlated inputs. A comparison of this approach 
with the older matrix-based formulation (Refs. 4, 5, and 6) showed 
the iterative method to be conceptually simpler, and probably 
more economical in computer time, especially as the number of 
correlated inputs increases. To verify this conclusion, it would 
be necessary to analyze an appropriate model by both methods then 
make a direct cost comparison. 

To find out whether the inputs to the system were correlated 
or not, it was necessary to calculate the OCF between the various 
inputs; if the value of the OCF was no greater than 0.1, then it 
was assumed that the two inputs involved were uncorrelated so that 
the source identification could be carried out in terms of OCFs. 

If OCF >0.1, the inputs were assumed to be correlated and the more 
complex PCFs had to be calculated. This test must be performed 
as a routine part of the analysis. It was found that, in some 
cases inputs designed to be independent were actually well corre- 
lated at certain frequencies due to structural feedback between the 
measurement points, ^^en the test indicates that the inputs are 
correlated at some (but not all) frequencies, the problem can be 
broken into two parts; a frequency basis and the appropriate 
coherence functions (ordinary or partial) calculated for each part . 
If the two frequency ranges are not well separated, a practical 
approach would be to assume that the inputs are correlated at all 
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frequencies and analyze accordingly. At frequencies where this 
is not true, the appropriate terms will drop out and the correct 
solution will be obtained, at the expense of some increase in com- 
puter time. 

The CF plots for this study used a linear scale for CF and a 
log scale for frequency, which is the format commonly found in the 
literature. For comparison, one plot was made with linear scales 
for both parameters, shown in Figure 32. Although the low fre- 
quency resolution suffers by this approach, the high frequency 
resolution is improved. For some applications, where the fre- 
quency range is fairly narrow, this format may be more desirable. 

The effect of increasing the number of data segments used for 
averaging was demonstrated, using data from Run 4. The input/ 
output OCF was plotted for 10, 50, 100, and 200 averages. It was 
concluded from variations in the plots that it is desirable to 
use 150 to 200 averages to ensure that the CF estimate is reason- 
ably accurate, without incurring excessive computer costs. Since 
the duration of the data segments affects the accuracy of the re- 
sults, as well as the number of segments, a trade-off between the 
two is necessary. A thorough investigation of this problem would 
be desirable, but was beyond the scope of the present study. 

Another factor affecting the accuracy of the CF estimates is 
the windowing process applied to the test data. In this study, 
only rectangular windows were used (i.e., the time history was 
terminated abruptly at the beginning and end of each record. Im- 
proved accuracy may be gained by using more sophisticated windowing 
techniques, but no optimum window exists for all applications and 
to find the best window for coherence functions would require 
further study. 


Phase II: Nonstationary Inputs 

The ordinary coherence functions between the nonstationary 
inputs and the panel response appeared to vary randomly as the 
time period spanned by the 20 segments used for averaging was 
increased. Even the first case, which was calculated over a data 
record of only 4 seconds, resulted in OCFs that were significant- 
ly different to the OCFs for the stationary part of the run. Also, 
the time variation of the inputs used for the test program was 
fairly slow, compared with the variation experienced during a 
spacecraft launch for example. In Reference 5, it is shown that 
spacecraft vibration during liftoff will typically double in RMS 
level in about 4 seconds. In the test program for this study, 
the RMS level doubled in about 20 seconds. 
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This indicates that to calculate functions from flight data 
with reasonable confidence, very high data sampling rates would 
be needed so that the time duration of the data records could be 
minimized without incurring excessive errors. Short-time averag- 
ing techniques discussed in Reference 5, could then be applied to 
calculate "locally-stationary" spectral density functions, and 
hence, a locally-s tat ionary version of the coherence function. 

The conclusions for the study, and recommendations for future 
studies, may be summarized as follows: 

Cl) For stationary data, it has been demonstrated that co- 
herence functions can be used to identify and quantify sources of 
dyanmic excitation, for both correlated and noncorrelated inputs. 

C2) The iterative approach formulated by Bendat for use on 
correlated-input systems is relatively simple to apply to a prac- 
tical system. 

(3) Bendat ' s approach is probably cheaper to run on a digi- 
tal computer than the earlier matrix techniques, but side-by-side 
comparisons need to be made. 

(4) Approximately 150 to 200 averages should be used in the 
computation of the auto- and cross-spectra used in the CF calcu- 
lation. 

(5) The use of windowing techniques on the measured data should 
be investigated, with the objective of reducing errors. 

(6) Even for relatively slowly -varying nonstationary data, 
the coherence functions cannot be calculated accurately by the 
usual stationary data analysis techniques. More work is needed 
to develop practical methods for handling nonstationary data. 
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